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Abstract 

We consider the nonlinear Dirac equation in one dimension (the massive Gross-Neveu model). We explicitly 
construct solitary wave solutions, and then study the linearization of the equation at a solitary wave. We present 
numerical simulations and justify them with explicit construction of some of the eigenfunctions. Then we present a 
WKB-based argument which justifies (but does not prove) the spectral stability of solitary waves of sufficiently small 
amplitude. We also compare our results with previously known numerical simulations. 

1 Introduction 

Field equations with nonlinearities of local type are natural candidates for developing tools which are then used for the 
analysis of systems of interacting equations. Equations with local nonlinearities have been appearing in the Quantum 
Field Theory perhaps starting with the articles by Schiff MSch51al ISch51bl . who studied the nonlinear Klein-Gordon 
equation in his research on the classical nonlinear meson theory of nuclear forces. These papers stimulated devel- 
opment of mathematical tools for tackling such nonlinear equations. Well-posedness for the nonlinear Klein-Gordon 
equation was addressed by Jorgens H Jor6 1 1 and by Segal [Seg63[. Nonlinear scattering was studied by Morawetz and 



Strauss in [Seg66 Str68 MS72|. Existence of solitary waves in nonlinear Schrodinger equation and Klein-Gordon 
equation was proved for a large class of nonlinearities in 0Str771 and 1BL83I . Stability of solitary waves takes its 
origin from IVK73I , where the (spectral) stability of the linearization at a solitary wave for a nonlinear Schrodinger 
equation has been studied. Orbital stability and instability of solitary waves in nonlinear Schrodinger, Klein-Gordon, 
and Korteweg - de Vries equations has been extensively studied in HSha831 |SS85l |Sha85l |Wei86l IGSS87I IBSS87I . 
The asymptotic stability of the standing waves in nonlinear Schrodinger equation was proved in certain cases in 
HWei85l IBP931 IS W92I ISW991 ICucO 1 1 IBS031 ICM08I . 

Much less is known for systems with Hamiltonians that are not sign-definite, such as the nonlinear Dirac equation, 



such as the Soler model BSol70l , and Dirac-Maxwell system. The existence of standing waves in the nonlinear Dirac 
equation was studied in [Sol70|, [CV86|, |Mer88|, and HES95I . The existence of standing waves in the Dirac-Maxwell 
system is proved in [EGS96| (for ui £ (— m, 0)) and [Abe98| (for uj £ (— m, m)). For an overview of these results, 
see flES02|. The local well-posedness of the Dirac-Maxwell system was considered in HBou961 . The local and global 
well-posedness of the Dirac equation was further addressed in HEV97I (semilinear Dirac equation in n ~ 3), BBouOOl 
(interacting Dirac and Klein-Gordon equations in n = 1), and in [; MNNO05 | (nonlinear Dirac equation in n = 
3). Numerical confirmation of spectral stability of solitary waves of small amplitude is contained in MChu07l . The 
asymptotic stability of small amplitude solitary waves in the external potential has been studied in IIBou08l . 



'Supported in part by the National Science Foundation under Grant DMS-0907968. 
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In this paper, we give numerical and analytical justifications (but not a rigorous proof) of spectral stability of small 
amplitude solitary wave solutions to the nonlinear Dirac equation in one dimension. 

ACKNOWLEDGMENTS The authors would like to thank Marina Chugunova for providing us with the numerical 
results on spectral properties of linearized Dirac equation (see also QChu071 ) which greatly stimulated our research. 
The authors are grateful to Dmitry Pelinovsky for many helpful discussions. 



2 Nonlinear Dirac equation 

We consider the nonlinear Dirac equation of the form 

idtip 



3=1 



x € 



(2.1) 



where 







N 



and 



with ip + = [-0!, . . . , ijjjf] being the Hermitian conjugate of ip. The Hermitian matrices a and (3 are chosen so that 
(— i Y?j=i a jdxj + (3m) 2 = (—A + to 2 )/, where A = Y^=i ®x * s me Laplace operator and / is the unit matrix. 
That is, a.j and (3 are to satisfy 



°4 = 7 > 



f3 2 =I, {a j ,a k } = 28 jk , {a„/3} = 0. 



We assume that the nonlinearity g is smooth and real-valued. We denote m = g(0). 

When n = 3, one can take Dirac spinors (N = 4 components). This equation with n = 3 and g(s) = 1 — s is the 
Soler model HSol701 , which has been receiving a lot of attention in theoretical physics in relation to classical models 
of elementary particles. The case n = 1 (when one can take spinors with N = 2 components) is known as the massive 
Gross-Neveu model HGN74I ILG751 . 

In terms of the Dirac 7-matrices, equation (12.1b takes the form 



iYd^ip = g(^tp)ip, where 7 = j3, -y J = 0ctj, d = d t , dj = d Xj 
The nonlinear Dirac equation can be written in the Hamiltonian form as ip = ~i5j,E, where 



-iip/3otjd X:j ip + ii), cHaji:) + G(^V) 



with G being the antiderivative of g such that G(0) = 0. Note that G{3£) = -F{SC) from 0SV86I . 
We consider the Dirac equation in K. 1 : 

i(d t + ad x )ip = g(iji/j)t3ip. 
As a and /3, one can take any two of the Pauli matrices. We choose 

a = -cr 2 , (3 = <x 3 , 

where the Pauli matrices are given by 



(2.2) 



dx, (2.3) 



(2.4) 



1 

°"1 = I 1 q I j Q"2 = I 

The nonlinear Dirac equation j2A\ takes the form 

i(d t + (-<T2)d x )ip = g(ip + a 3 ilj)(T 3 t(j 

We can rewrite this as a system 



— i 

1 



1 

-1 



idtip! = d x tp 2 + g(\ipi 



id t ip 2 = -d x ^ l -g{m 2 -m 2 )^ 2 



(2.5) 

(2.6) 
(2.7) 
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3 Solitary wave solutions 

We start by demonstrating the existence of solitary wave solution and exploring their properties. 
Definition 3.1. The solitary waves are solutions to (12. U of the form 

S = {t/>(x,t) = Mx)e~ ZUJt ; K G H\R n ,C N ), u G R}. 
The following result follows from ICV86I . 
Lemma 3.2. Assume that 

m = g(0) > 0. 



(3.1) 



Let G be the antiderivative of g such that G(0) = 0. Assume that for given ueR,0<u<ra, there exists X u > 
such that 

cjX^ = G(X U ), to ^ g(X u ). and ujX < G(X) for X G (0, X u ). (3.2) 
Then there is a solitary wave solution tp(x, t) = <p UJ (x)e~ lu ' t , where 



v{x) 
u(x) 



v, u G H 1 ^), 



(3.3) 



with v, u real-valued. 

More precisely, let us define X(x) and W(x) by 



X = v 2 - u\ 
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Then X(x) is the solution to 

X" = -d&(-2G(X) 2 + 2cj 2 X 2 ), X(0) = X u 
andW{x) = -±X'(x). 
Proof. From ( 12.71 ). we obtain: 



X'(0) = 0, 



(3.4) 



(3.5) 



ojv — d x u + g(\v\ 2 — M 2 )i>, 
lou = —d x v — g(\v\ 2 — |u| 2 )u. 

Taking into account that both v and u are real-valued, we can rewrite ( 13.6b as the following Hamiltonian system, with 
x playing the role of time: 

— d x u = —cjv + g(v 2 — u 2 )v = d v h{4>), 
d x v = —lou — g(v 2 — u 2 )u = d u h(4>), 



(3.6) 



(3.7) 



where 



h{cp) = -^{v 2 + u 2 ) + \ G {v 2 -u 2 ). 



(3.8) 



The solitary wave corresponds to a trajectory of this Hamiltonian system such that rim :r _ ( ± o v(x) = lim x ^± oc u(x) = 
0. Since G(X) satisfies G(0) = 0, 

h(v(x),u(x)) = 0. (3.9) 

Thus, 

iu(v 2 +u 2 ) = G(v 2 - u 2 ). (3.10) 
Studying the level curves which solve this equation is most convenient in the coordinates 

X = v 2 -u 2 , 2f = v 2 +u 2 ; 

see Figure[3] We conclude from (13.10b and Figure[3]that solitary waves may correspond to < m, u> ^ 0. If u> > 0, 
then the solitary waves correspond to nonzero v, with u changing its sign. 

Remark 3.3. If u> < 0, then the solitary waves correspond to u ^ 0, and to v changing the sign. 
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Figure 1: Existence of solitary waves in the coordinates 3t~ = v 2 — u 2 , 3? = v 2 + u 2 . Solitons with lu > and uj < 
correspond to the bump on the v 2 axis and to the dotted bump on the u 2 axis (respectively) in the first quadrant. 

The functions S'(x) and & (x) introduced in ( 13.4b should solve 

w = -(v 2 + u 2 ) g {^) + uj,r = -±G(&)g(ar) + Lu$r, ^ ' 

with the limit behavior limui^oo 3£(x) = 0, lim^i^^ W{x) — 0. In the second equation in (13.1 IK we used the 
relation (13 - Sb and the identity h(4>) = 0. The system ( 13.1 It can be written as the equation on 3£: 

3C" = _a 3r (_2G(JT) 2 + 2w 2 Jf 2 ). (3.12) 

This equation describes a particle in the potential U UJ (^') = —2G(3£) 2 + 2uj 2 3£ 2 \ see Figure|2] Due to the energy 
conservation (with x playing the role of time), we get: 

aril ar/2 

^--2G( < r) 2 + 2c J 2 jr 2 = ^— + u u (ar) = o. (3.13) 

Using the expression for 3£' from d3.1 lb . relation (13.13b could be rewritten as 

= + U^{3C) = 8oj 2 ^ 2 - 2G(,T) 2 + 2uj 2 SC 2 = 2uj 2 {Av 2 u 2 + {v 2 - u 2 ) 2 ) - 2G 2 , (3.14) 

which is equivalent to (13.10b . 

For a particular value of uj, there will be a positive solution 3£{x) such that lim x ^±oo 3j(x) = if there exists 
> so that ( 13.2b is satisfied (see Figure|2]i. For convenience, we assume that JT(0) = 5E W , so that S£{x) is an 
even function. 



U U (3T) = -2G 2 (,T) + 2oj 2 % 2 

= — 1 ^^***^^cj 



Figure 2: Effective potential U u (^). A solitary wave corresponds to a trajectory which satisfies JT(0) = SK^, 
lim| x |_ 0O 3£(x) = 0. 

After we find 3£(x), the function & (x) is obtained from the first equation in ( 13.1 lb . Knowing 3£(x) and & (x), 
one can express v(x) and u{x). □ 
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4 Explicit solitary waves in a particular case 

As shown in 1LG75I for the massive Gross-Neveu model, in the special case of the potential 

Of 2 

G(X) = X- — , 



(4.1) 



the solitary waves can be found explicitly. Substituting G(X) from (14.11 ) into ( 13.13b . we get the following relation: 

dX dX 



dx 



We use the substitution 



' 2^/(X - ^" 2 /2) 2 -uj 2 X 2 2&y/(l- X/2) 2 -oj 2 
X v 



cos 29 ' 



X = 2 1 



Then 



dx 



dX 



) 2m sin 26 

' cos 2 26 



cos 26 



dS 



dO 



2XJ(1-X/2) 2 -uj 2 am uj \ r~Tp~ 

VV ' ' 4 U ^2e>y cos 2 26 

y 7 /! + tan 9 



cos 28 — uj 



x = — In 
2k 



where 



Also note that 



and then 



K=y/T 



UJ 2 , 



yJJL — tan 8 ' 
M ~~ 1+uj' 



(y / ]2 + tanO)e 2KX = ^JJl — tanO, tan© = — y^tanhra. 



X = 2 1 - 



i 1 cos 29 ) 



= 2 1- 



2 cos 2 6 - 1 



= 2 1 -cj 



1 + tan 2 6 
1 - tan 2 6 
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4w 



-X' 



i * d6 

4 cos 2 26 1 ' dx 

X 



1 2 



Denote 
Then 



2T(a:) 



cos 28 

There are the following relations: 



tan 26 - 
2 2 1 - tan 2 6 



X[x) = v 2 (x)+u 2 (x). 

1 + tan 2 6 



(-2sin26)(cos26-cj) 
X- 



4 cos 2 26 
■SBT 2 tan 6 y/TItanliKa; 



cos 26 



1 - tan 2 6 



1 -uj 



1 — /i tanh kx 



1 + tan 2 6 
1 - tan 2 6 



v = s/W cos 6, u = —yfX sin 6, 

1 



Jf = J"cos2e, <3f 



-2? sin 26. 



(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 
(4.7) 

(4.8) 

(4.9) 
(4.10) 

(4.11) 



(4.12) 
(4.13) 



Remark 4.1. By ( 14. 1\ , tan 6 changes from JJi to — ^JJl as x changes from — oo to +oo. Thus, in the limit u> — > 1, 
when /x -> 0, one has X w JT, while |^| < 3?\//Z. 
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5 Linearization at a solitary wave 

To analyze the stability of solitary waves we consider the solution in the form of the ansatz 

u(x,t) = (4> u (x) +K{x,t))e- luJt , <j>„ 
Then, by flTfl . 

iK + ujK = -i^2ad Xj K + (3 [g{{<f>«, +H){<l>u +n)){<j> u +11) - s(<^<M<?M 

3 

The linearized equation on 7Z is: 

ill = -iJ2 ad ^ n -uK + (3 [g$MK + $>JR. + R^yftWi] . 



(5.1) 



Let us write this explicitly. Namely, let 1Z = 
takes the following form: 



Pi + Kl 
P2 + K2 



, with pj, <Tj real-valued. The linearized equation on p, c 



dtp = d t 
where J corresponds to 
We can write (15.2b as follows: 



Pi 






' P 2 




Pi 




Pi 


P2 


"'I 




-Pi 


— o; 


P2 


+ .9 


-P2 


Si 












?1 


^2 






-4 




? 2 







+ 2(^/91 - Wp2)ff' 



w 




— u 







j 








p = JJlfuP 



h 

-h o 



J 



g - lo 3j. 
-5 X -.g - w 



J 2 
-J 2 o 



.9 - w d x 
—d x -g - lo 



2g' 



v 

U 



uv 

2 








where I is the 2x2 unit matrix. 
We define M± as follows: 



g - W 9as 

-9,, -o - 



yc + 



2g'v 2 + g — lo d x — 2g'vu 
—d x — 2g'vu 2g'u 2 — g — lo 



Then the linearization at the solitary wave takes the form 

dt 



p 




I 2 ' 











p 











p 


<r 




I 2 

























Remark 5.1. A similar block-decomposition was obtained in |CP061IChu07l| . 

Remark 5.2. Alternatively, to analyze the stability of solitary waves, we need to consider the operator 

& = (</>„) -wQ"(</> u ). 
Let <P = (^1,^2), where <Pi = tpi + $2 = <p 2 + 1X2, where cpj and Xj are real-valued. We have: 



Qm = \ I (fi + + xi + X2) Q"(*) = U 

1 Jm. 



(5.2) 



P, (5.3) 
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The energy functional is given by 

= / (-«'#* a<2>' + G($$)) dx 



1 
-1 



1*2 



dx 



({{<Pi ~ «Xi)(V2 + *X 2 ) - (V2 - «X2)(vi + iXi)) + G(<Pi + xl- ( fil~- xD) dx. 
For its second derivative, we have: 

fiXi 






9z 








■d x 




















0, 








-d x 






Evaluating this at the solitary wave, # 
Then 

5.1 Spectra of J?± 



2G" 



-<Pl<P2 
flXl 



-VlX2 
V2X2 



V2X1 

Xi -X1X2 

V1X2 V2X2 -X1X2 X2 



■^2X1 
P2X2 

(v, u), with both v, u real, we have: 



G' 






9a 




























9a 








-9a 






v 2 










-vu 



































2G" 



,yf = E"{<j>) - wQ"(<£) = £"0) - u)I 4 



+ G' 



Jt?+ 
JffL 





-1 

1 







1 
1 




While we are ultimately interested in the spectrum of the operator L u , we start by analyzing the specta of J$?± which 
are easier to compute and which will shed some light on the behaviour of the full operator . 
Define 

m± = 1 ± lu. 



Lemma 5.3. <r cont (Ji?±) = R\(— m+, m_). 

Proof. The symbol of at x — ► ±00 and £ = is given by lim|.j.|_ ! . 00 (j{J%±){x, £)L =0 
eigenvalues correspond to the edges of the continuous spectrum. 

Lemma 5.4. a d {,3^-) D {-2w,0}. 



771_ 

— m + 



Its 
□ 



Proof. By (ED, J#i 



= 0, hence = 



£ ker . Moreover, there is a relation 



= 0, 



u 




g + ui d x 




u 


V 




-d x -g + oj _ 




V 



where again we used ( 13.6b . 

The numerical computations show that there are no other eigenvalues in J$?_ ; this agrees with IIChu07l . 
Lemma 5.5. <7 d {M+) D {-2w,0}. 
Proof. One can immediately check that 



(5.4) 
□ 



Moreover, 
The last equality is due to 







u 




V 





' v' ' 









2g'v 2 + g — uj d x — 2g'vu 
—d x — 2g'vu 2g'u 2 — g — uj 



v 
u' 



2g'v 2 + g + uj d x — 2g'vu 
—d x — 2g'vu 2g'u 2 — g + uj 



ft 

V 



g + uj d x 
-9a -g + lu 





u 




V 



□ 
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Figure 3 : Spectra of H+ (left) and of H — (right) versus u>, computed in IChu07l via Scalapack library on a parallel computer cluster. These 
plots are courtesy of IChu07l . 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 

Figure 4: Our recalculation of a(H+) via the shooting method. Note the presence of the eigenvalue A = 2lu of cr(H+); the eigenfunction 
corresponding to A = 2lo has been found in the explicit form, which seems to confirm our plot. Bifurcation at u> 0.5 is absent, which could be 
due to a better resolution of the shooting method. We do not know the reason for the discrepancy in computations of cr(H-\-). We would also like 
to point out that our recalculation of cr(H— ) is in perfect agreement with the right plot on Figuref5]obtained by .('lui()7|. 



The rest of the spectrum of Jt? + is computed numerically; see the spectrum of the operator H+ ;= on 
Figure HI (We consider — Jtif + for easier comparison with the spectra of H + and L u from BChu071 represented on 
Figure [3]) 



5.2 Spectrum of L u 
Lemma 5.6. 1. 

2. 

3. a d (L) 3 {±2iu}. 



Ocont{Lui) = (—00, -1 + Lj] D [1 - CJ, +00). 

ker = Span • 



d x <j) 




' " 










Proof. Let tp 



ii 

v 



Due to the above results, 



±iip 



= -2lu 



Hip 



To find the continuous spectrum of L u , we consider the limit of L u as x — > ±oo, substituting v, u by zeros, g by 
g(0) = m = 1, and d x by ££. We would also like to find the asymptotics of functions which solve (L — A) "J = for 



a particular A G C. We assume that such a function has the form ^ 

(L-X) 



At; 



±00, we have: 



r 




-A 


H 




r 


s 




-H 


-A 




s 



with H = 
where H 2 = 



—it; — m. 
m 2 _ + £ 2 -2iuj£ 
2iuj£ ml + £ 2 



coinciding with both Jt?± at \x\ — ► oo. Hence ffs = Xr, Hr = —As, iJ 2 r = — A 2 r, 
Therefore, A satisfies 



= det(i? 2 + A 2 ) = (ml + £ 2 + A 2 )(m 2 + £ 2 + A 2 ) - 4w 2 £ 2 , 
(1 + lu 2 - 2uj + £ 2 + A 2 )(l + uj 2 + 2lu + + A 2 ) - Aoj 2 i 2 = 
(1 + oj 2 + C 2 + A 2 ) 2 - 4w 2 - 4w 2 C 2 = 0, 1 + £ 2 + w 2 + A 2 = ±2uVl + f 2 , 



1 + f + uj 2 ± 2w V 1 + C 2 = -A 2 , w± v / l + C 2 = ±iA 



Therefore, 



A = ±i(w± v^l + £ 2 ), ±iA = w± ^1 + ^ 2 , V 7 ! + C 2 = ±w ± iX, 
The continuous spectrum of L u (values of A which correspond to £ £ M.) is 

VcontiLu) = (—too, — «m_] n [i77i_, -Moo). 



e±± =±v / (^±*A) 2 -l. 



□ 



Possible resonance points (£ = 0): 



A 4 + 2A 2 (1+w 2 ) + (cj 2 -1) 2 = 0, 



A 2 = -(1 + LU 2 ) ± y/{l+L0 2 ) - (1-LJ 2 ) 2 = -l-L0 2 ± 2UJ, 

A = ±i(l±cu). 

For the eigenfunctions of L w corresponding to a particular A € C one can take 

2iw£ 
mi + £ 2 + A 2 



r = 


n 


= A 




r 2 





(5.5) 



which is an eigenvector of H 2 corresponding to —A 2 . 

Remark 5.7. This definition of r gives an eigenvector of H 2 corresponding to — A 2 everywhere except at A = ±im_ 
where (15.5b vanishes. Then one could take 



r = 


ri 


= A 




r 2 





m\ + £ 2 + A 2 
-2iaj£ 



(5.6) 



which represents an eigenvector of H 2 corresponding to the eigenvalue — A 2 everywhere but at A = ±zm+ (where 



vanishes). 

Then we use the relation Hr = — As to find 



S\ 












S2 






— m + 




m 2 _ 



2iu$, 

_1_ tf2 , \2 



i£(2m-uj + m 2 _ + ( 2 + A 2 ) 
2w£ 2 - m+(m 2 _ + £ 2 + A 2 ) 



Denote 



C±± = ±£± = ±\/(w±a) 2 -i = ±v/(-aTw) 2 -i. 

Let A = a + ife. We assume that both a and & are non-negative. Then 



f±± = ± V / ((6-ia)Tw) 2 -l = ±V(6 T ^) 2 - 2ai(& T w) - a 2 - 1. 
For s/z, we choose the branch of the square root such that Rc s/z > 0. Define 

C±± = {zeC; ±Rcz > 0, ±Jmz > 0}. 



We expect that the eigenvalues emerge from i(l + lu). Assume that 1 < 6 < 1 + a;, 0< a< 1. Then 



£+_ = + v / {b + lu) 2 - 2ai(b + uj) - a 2 - 1 « +^(1 + 2lo) 2 - 2az(l + 2lu) - 1 G C+_, 



£ ++ = + \/(& - w) 2 - 2ai(b - u) - a 2 - 1 ~ + \/-a 2 - 2ai = —iy/ 'a 2 + 2ai £ C+_, 

-+ 



e-+ = g c_+, e+ = e++ g c__. 



The eigenvector corresponding to £ = £+_+_ is given by 



n 




^2 




Si 




_ s 2 





\2iuj£; 

\{m 2 _+e + \ 2 ) 
-2m-iu(, - ii{m 2 _ + £ 2 + A 2 ) 
-2lu£ 2 + m. + {m 2 _ + £ 2 + A 2 ) 



Noting that 



we have: 



mi + g ± + A = (1 — lu) z + (lu± iXY — 1 + A = 2lu - 2uj ± 2i\u = 2u){-m- ± 2iA), 



iA — m_ 
A + im_ 



it 

—i\ — m_ 

-£ 

A — !TB_ 



5 — - + 



iA — m_ 

HE 

A + im_ 



— ?'A — m_ 

A — !TB_ 



£s in the upper half-plane correspond to waves decaying as x — » +oo 









• £++ • £+- 



£s in the lower half-plane correspond to waves decaying as x — > — oo 

Figure 5: Im^ + ± < 0, hence e l ^ +±x decays for x — > — oo. On the other hand Im£_± > 0, hence e^- ±a: decays for 
x — > +oo. 

Denote by X~ the subspace in C(M, C 4 ) with even first and third components and with odd second and fourth 
components. Similarly, denote by X + the subspace in C(R, C 4 ) with odd first and third components and with even 
second and fourth components. Then C(M,C 4 ) = X + © X~ . Noticing that L u is invariant in X + and X~ , we 
conclude that all eigenvalues of L u always have a corresponding eigenfunction either in X + or in X~ (or both). To 
find eigenvalues of L u corresponding to functions from X + , we proceed as follows: 

• For A € C, construct solutions 1 < j < 4, to the equation L^S = X^> with the following initial data at 

x = 0: 



Then * 3 g X~, while * 2 , v[/ 4 g 
• Pick vectors S |_ and S such that (L a 



A)H_±e^- ±x = 0, where Im£_ ± > 0. 



Define the analog of the Evans function E (A, R) = dot 



Similarly, define E + (A, i?) = det ^ 2 L =fl , *4L =R , 5 



10 



• The condition E (A, R) — means that a certain linear combination of and ^3 G X~~ (or, respectively, ^2 
and ^4 G X + ) coincides at x = R with a decaying solution of [L x — +00 — A)\& = 0. Therefore, if A G crd(L), 
then either liiriR-^+oo -E" (A, i?) = or IiraR^+oo R) = (or both). 

Numeric computations do not show any unstable eigenvalues between to = 0.2 and ui = 0.9 (we do see exact 
eigenvalues at A = and A = 2uii). On Figure|7l we plot the results of our computations for oj = 0.2, in the vicinity 
of the point where - presumably - an unstable eigenvalue from IIChu07l (see Figure|6l right) was to be located. 



6 Resonances of Jif + 

We consider the operator 



2g'v 2 + g — lu d x — 2g'vu 
—d x — 2g'vu 2g'u 2 — g — ui 



6.1 Resonance of at A = m_ = m — u 

The resonance eigenfunction (M, N) is to satisfy the equation 



" —N' ' 




' -3v 2 + u 2 


2uv 




' M 


M' 




2uv 


-2 + v 2 - 3u 2 




N 



We are looking for a solution of the form 

lim 



x — > — oc 



M{x) 
N(x) 



lim 

x — - + oc 



M(x) 
N{x) 



±1 




Assuming that u> < 1 and keeping the principal contribution (see Remark l4~TT i, we get: 

-N' = -3v 2 M, M' = (-2 + v 2 )N, 
M" + (2- v 2 ){3v 2 )M = 0. 

The WKB approximation for a solution is: 



M{x) « cos ( f ^(2-v 2 {y))3v 2 {y) dyj . 



The resonance condition takes the form 



6v 2 dx w Vo / \J 3£ [x) dx = nn, n G 



Since 3£ m v 2 m 2(m — uj)e 2k N, where k = \/m 2 — uj 2 , we have 



(6.1) 



V6 / VWdx w \/6 / dxfnV6 ^2(m~ uj)er Klxl dx w \/6 ■ 2 



^2(m - w) _ 4V3 



y/m 2 — lj 2 y/l + u 

Remark 6.1. Let us check the applicability of the WKB method. As a matter of fact, the applicability condition 
3£(x) << 2£(x) does not hold. Let us check the magnitude of the error. We consider the equation 

2 



-Z" 



z = o, 



L 



cosh x > 

which has the exact solution Z(x) = — tanhx. The WKB method gives the approximate solution 

'■ x V2 



J WKB 



(x) = cos ^ 



coshy 



dy 
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Figure 6: Spectra of Luj at oj = 0.2 (left) and to = 0.7 (right), computed in fChu07j via Scalapack library on a parallel computer cluster. Note 
the absence of exact eigenvalues A = ±0.42 on the left plot. These plots are courtesy of | Chu07 '. 
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Figure 7: Trying in vain to corroborate an unstable eigenvalue of Luj at oj = 0.2 by computing the Evans functions ^^(A, R) (with R = 20) 
for even-odd-even-odd spinors from X~ (left) and odd- even-odd- even spinors from X + (right) in the vicinity of A ~ 0.15 + 0.82. The unstable 
eigenvalue III (A pa 0.15 + 0.82) from Figure [6] (left) would be somewhere in the middle of the above plots. The regions where the Evans function 
takes positive real values are marked with "x", while positive imaginary values are marked with "o", An eigenvalue at a zero of the Evans function 
would be at the intersection of borders of the marked regions. None of such intersections takes place on either plot. 
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1. The total phase change for tanhx as x changes from — oo to +00 is tt, while total phase change of Zwkb is 
L dy = V2n. 

2. The asymptotics for x < 0, \x\ 3> 1 for the solution = — tanhx and its approximation Z^/kb^) are 

— e~ x 1 — e 2x 



= — tanhx = 
Zwkb(x) w cos 



+ e 



1 + e 



2.r 



1 - 2e 



2.r 



2.r 



„ ) = cos ~ 1 — . 

2 V V^2/ 4 

We see that the WKB asymptotics are imprecise. 

The numerical simulations show that the "phase shift" by 2tt (n = 2 in (16. jj ) corresponds to the limit w — ► 1. At 
the same time, the phase shift which ( 16. Il l gives for w = 1 is approximately 4-^/3/2 ~ 4.9. Numerical computations 
show that the phase shift by 37r (n = 3) corresponds to lu w 0.367, the phase shift by 4-7T (n = 4) corresponds to 
w f» 0.205. For these values of u>, our approximate condition ( 16.1b becomes very imprecise since it was derived for 

<*> < 1. 

6.2 Resonance of Jff + at A = — m + = —m — uj 

The resonance eigenfunction (M, iV) is to satisfy the equation 



" -N' ' 




M' 





We are looking for a solution of the form 



lim 

x — > — 00 



" M(x) ' 




' " 


N(x) 




1 



2 - 3i> 2 + u 2 2uv 
2uv v 2 — 3u 2 



lim 

x — >+oo 



M 
N 



' M(x) ' 







N(x) 




±1 



Keeping the principal contribution, we get: 

—N' = (2 - 3v 2 )M, M' = v 2 N, 
N" + (2 - 3v 2 )(v 2 )N = 0. 

The WKB approximation for a solution is: 

N(x) « cos ( / ^/(2~3v 2 (y))v 2 (y)dy) . 
J —00 

The resonance condition takes the form 

/ dx^V2 \J 3£{x) dx = nn, n e Z. 
Js. Jm 

7 Resonances of L w at A = im± 

Recall that m± = m ± u>. The eigenfunctions corresponding to A = i(m ± ui) = im± satisfy 

im± 



p 







h ' 








P 


c 




I 2 














(6.2) 



im±p = Jtf'-q, im±q = —Jtf+p. 
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Explicitly, 

im±p = JtfLc, 



g-uj d x 
—d x —g — uj 



2g'v 2 + g — oj d x — 2g'vu 
—d x — 2g'vu 2g'u 2 — g — lo 



p'l = —2g'vupi + (2g'u 2 — g — w)p2 + «m±ft, 
— p' 2 = {2g'v 2 + g - u)pi - 2g'vup 2 + im±?i, 
?1 = (-.9 - w)ft - im±p 2 , 
-?2 = (.9 ~ w Ki - im±p 1 . 



Keeping the main contribution, we get: 



p'l = (~m + + ££)p 2 + im±^ 2 , 
—p' 2 = (m_ - Z3?)pi + im±gi, 

= (— m + + <^")ft - im±p 2 , 
-?2 = (m_ - - im±p 1 . 



7.1 Resonance of L„, at A 



We take: p 2 , ft, ft) 



pi = (— m+ + ^")p2 + im + q 2 , 
—p' 2 = (m_ — 3^T)/>i + im+ft, 
^ = (— m+ + <3T)ft - im + p 2 , 
-q 2 = (m_ - - im+p!. 

(0, 1, 0, — i). We differentiate the second and fourth equations and disregard 3£'(x): 



-p 2 = (m_ — 3^)^ + im + <^ = (m_ — 3^")((— m + + 3£)p2 + im+ft) + + ^)ft — i m +P2) 

-q 2 = (m_ - - im + p[ = (m_ - JQ((-m + + $T)ft - iTn + p 2 ) - im + ((—m + + 2£)p 2 + im+ft), 

which takes the form 

p' 2 ' + rn^_/52 + (m,- — 3^')(— m + + %}p2 + (m- — 3^)im+^2 + im+(— m+ + =5£~)ft = 



?2 ^ ("i- — J%~)im+p2 — im + (—m + + 3fc)p2 + ™+ft + — m+ + <5T)ft = 

Disregarding terms with S^ 2 , we get: 



Pi 


" + { 


ft 





+ < (mA — m_m_|_) 



1 -i 
i 1 



3m + + m_ — 2im+ 
m+ + m_ 



} 


P2 




ft 



= 0. 



Denote T2 = /?2 + *ft, T 2 = P2 — K2- Coupling the above system with (1, i), one gets 



Pi 
ft 



r 2 + 3K"[2>m + + m_, —im + + im_ 
Coupling the above system with (1, — i), one gets 

T 2 +2m_)_(m-)_— m_)T2+=^"[3m++m_, —2>i r m + —im 



= r' 2 ' + ,3f;{m+ + m_)r 2 + 2 Jfm + T 2 = 0. 



(7.1) 



P2 

ft 



= T2+2m+(m+— m_)T2+>^"(3m++m_)T2 = 0. 



Thus, T2 = 0. The WKB approximation of ( 17. U yields r 2 (x) = C cos(J^ \j2m,%\y) dy). 

Since the solution (p, <r) can be broken into odd and even components, the resonance condition takes the form 



\/2ml"(a;) dx 



n € Z. 



(7.2) 



Remark 7.1. Let us note that the resonance condition ( 16.2b for ^+ at A = — m+ coincides with the resonance 
condition ( 17.2b for L w at A = ±1771+ (both derived under the same approximations), suggesting similar bifurcation 
pattern for at — m+ and for L w at ±zm+. Although these criteria in present form are very imprecise, the absence 
of bifurcations in the spectrum of J4° + indeed matches the absence of bifurcations in the spectrum of L u . 
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